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Abstract 

Many processes in nature seem to be entirely controlled by transition rates and the corresponding statistical dynamics. Some 
of them are in essence quantum, like the decay of excited states, the tunneling through barriers or the decay of unstable 
nuclei. Thus, starting from first principles, those systems should be analyzed on the basis of the Schrodinger equation. In 
the present paper we consider a two level system coupled to an environment which is basically described by an two-band 
energy scheme. For appropriately tuned environment parameters, the excitation probability of the two level system exhibits 
statistical dynamics, while the full system follows the coherent, unitary pure state evolution generated by the Schrodinger 
equation. 
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1. Introduction 

Quite naturally the question for the decay behavior 
of quantum states is not new and thus there is a con- 
siderable amount of theories aiming in that direction. 

First of all Fermi's Golden Rule should be men- 
tioned, a formula of enormous practical importance. 
Nevertheless its derivation is based on the assumption 
of a short perturbation and thus in can hardly account 
for a full, continuous decay process, all the way down 
to equilibrium. 

The decay process associated with the spontaneous 
emission of an atom is described by the Weisskopf 
Wigner Theory in a way that is not based on a short 
perturbation. But this theory describes particularly an 



Corresponding author. 

Email address: jgemmer9uos.de (Jochen Gemmer). 



atom coupled to the electromagnetic field in the vac- 
uum state and is thus hard to generalize to arbitrary 
environments in arbitrary states. 

In the context of the system-environment scenario 
there are also the theories based on quantum master 
equations. Most of them involve some projection op- 
erator technique as well as the Born approximation 
and typically systems must be Markovian [1]. Many 
derivations are based on an initial state which is a 
product state with a thermal bath part. Some deriva- 
tions even assume this structure for the full relaxation 
process [2,3]. This assumption is often backed up by 
the argument that since for no interaction no system- 
environment correlations could be generated, weak in- 
teractions (Markovicity) could only generate negligi- 
ble correlations [4,5]. Furthermore, since those master 
equations result in maps on the considered system of 
the form p(t + r) — M.(r)p(t) that allow for an itera- 
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tion as M(ti + T2) = M(ti)M(t2) [6], the final state 
with regard to one map might be the initial state with 
regard to another. This fact also seems to support the 
idea that at least weak interactions could not produce 
considerable correlations. 

In the following we thus shortly comment on the 
question of relaxation and correlations in general, and 
then introduce an alternative theory of the relaxation 
process that takes arbitrary pure, possibly correlated 
or even entangled full system states into account. 



2. Relaxation and Correlations 

Is it possible that a full bi-partite system (system- 
reservoir-model) undergoes a unitary transformation, 
such that the purity of the considered system decreases 
(entropy increases, relaxation), without substantial 
system-reservoir correlations being generated? 

This question shall be addressed within this section. 
To those ends we specify the "correlations" , p c as an 
addend of the full system density matrix, p: 

p c := p-p s ®p r (p 3 := Tr r {p} , p r := Tr s {p}) (1) 

Obviously p a (g) p r specifies the uncorrelated product 
part of the density matrix. To measure the "size" of 
the correlated and the uncorrelated parts we use the 
absolute value, P, of an operator: 



Px ■= \/Tr{p2} :r = s,r,none (2) 

Evidently P is also the purity of the corresponding 
(sub)system. To decide whether or not correlations 
are negligible we finally want to consider the "correla- 
tions/product" coefficient rj: 

Pc 



V 



Ps Pr 



(3) 



If r\ <C 1, correlations may safely be neglected. 
Computing the size of the correlations yields: 

P 2 =p 2_ m {/3/5s/3r} + P 2p2 (4) 

Since the trace of a product of two hermitian matrices 
fulfils the conditions on an inner product, one finds 
through application of Schwartz' s inequality: 



\Tr{ P Pspr}\ < PPsPr 

Inserting this into (4) yields 

Pc > (P-PsPrf 



(5) 



(6) 



or for the coefficient 77 



P 

TlPr 



(7) 



P is invariant under unitary transformation. Often the 
reservoir is assumed to be exactly stationary which 
might not precisely hold true, nevertheless P r (0) ~ 
P r (t) should be a reasonable approximation for large 
reservoirs. Thus the only quantity that may substan- 
tially change upon relaxation on the right hand side of 
(7) is P s . And if P s decreases, r\ obviously increases. 
For the case of a stationary bath and an initial product 
state one finds 

Ps(0) 



Ps(t) 



- 1 



(8) 



This lower bound on r\ may easily take on rather high 
values, e.g., for an N- level system coupled to a bath 
in the high temperature limit (kT much larger than 
the level spacing) one gets for an initially pure state, 
Ps(0) = 1, 

r)>VN-l (9) 

This result is absolutely independent of the interaction 
strength. It only connects a decrease of purity (increase 
of entropy) to an increase of system-reservoir correla- 
tions, regardless of the timescale on which this relax- 
ation process happens. Thus we conclude that, quite 
contrary to the idea of system and bath remaining un- 
correlated correlations are generically generated upon 
relaxation. 



3. The Hilbert Space Average Method 

The method of describing the relaxation process we 
are going to present here relies on a step- wise solution 
of the Schrodinger equation. Naturally, this approach, 
like all other approaches, relies on approximations. 

The crucial approximation here is the replacement 
of some specific quantities by their Hilbert space aver- 
ages. Hilbert space averages are averages of quantities 
defined as functions of the pure full system state, over 
sets of such states which share a crucial common fea- 
ture, like, e.g., in our case, the same excitation proba- 
bility for the considered system. If the distribution of 
excitation probabilities among the states of the respec- 
tive set is broad, the replacement of an actual value 
by its Hilbert space average will most likely be a bad 
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approximation. If the distribution is tightly centered 
around its mean value, the replacement still only rep- 
resents a "best unbiased guess" but with distributions 
getting narrower, this guess obviously gets better and 
better. In our case there is evidence that the width of 
the distribution vanishes with an increasing number of 
environment states. 

Based on this idea we derive a rate equation, that 
describes the decay which is controlled only by the 
Schrodinger equation, up to some fluctuations. 



4. The System-Environment Model 

The energy scheme of the situation we are going to 
analyze is depicted in Fig. 1. A two level system, g 
("gas"), is in contact with a "many level" environment 
or "container" , c (We use this nomenclature for purely 
historical reasons [7]. Only the relevant parts of the 
spectrum of the environment enter the model. These 
are, in this case, two "bands" of width 5e, containing ATf 
(No) equidistant eigenstates in the upper (lower) band. 
Therefore the level spacing within the upper (lower) 
energy "band" is 



"gas" g 



"container" c 



1 Nf 



and 



(10) 



In the following, quantities of the "upper band" of the 
environment get the subscript 1, where as quantities of 
the "lower band" get the subscript 0. If we, e.g., con- 
sider an evolution from an initial state, with the sys- 
tem in the excited state 1 1) and the environment in the 
"lower band" , 0, due to overall energy conservation the 
only other set of states that the full system can evolve 
into, is the set with the considered system in the ground 
state |0) and the environment in its "upper band", 1. 
The Hamiltonian within the relevant subspace of the 
entire Hilbert space may thus be organized as follows, 
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Fig. 1. Discrete two-level system coupled to a 
quasi-continuous container system. This set-up exhibits, for 
a sufficiently high state density in the container system, 
and an adequately tuned coupling, an exponential decay of 
an excitation in the gas system. 

where i (j) count the levels in the upper (lower) "band" 
of the environment. The Hamiltonian is displayed in 
the eigenbasis of the uncoupled system, for simplic- 
ity we assume for the moment that the coupling V 
only adds terms to the Hamiltonian in the off-diagonal 
blocks. This corresponds to an coupling which may give 
raise to energy exchange between the subsystems. 

We now introduce two projectors, which project out 
the upper (lower) part of the state of the system 

P ox := ® l (c) , P gr :=|0)(0j®i (c) , (12) 

where i^ c ' is the 1-operator in the environmental sys- 
tem. In the following we call that part of the wave vec- 
tor that corresponds to the considered system in the 
excited state \ip e *) and the part that corresponds to 
the system in the ground state \tp &T ), i.e., 

|^ cx ) :=P™\4>), \4>n :=P gl M (13) 

=> |V) = |^ CX ) + |^ Sr ). 

Note that neither |i/>° x ) nor \i> gI ) are normalized indi- 
vidually. 

To analyze this model we first transform to the Dirac 
or interaction picture [8] 

U(t,0) := U := e-t" ot , |V>i) := U$\1>) , (14) 
Vj :=U%VUo, 

where Ho is the Hamiltonian of the uncoupled system. 
The Schrodinger equation in this representation reads 



(15) 



where both states and operators are now time-depend- 
ent, i.e. also V\ is a time dependent operator, but pre- 
serves the off-diagonal block form as before. 

The crucial quantities in the context of a decay to 
equilibrium are the probabilities to find the system in 



3 



its exited (ground) state, W eyi (W gr ). Due to the diag- 
onality of Ho those quantities have the same represen- 
tation in the interaction as well as in the Schrodinger 
picture, 

vr x = (VTIVO = W> cx lv> cx ) (i6) 

w sr = <v> g w> = {^\4>n ■ 

For simplicity we omit in the following the interaction 
picture subscript "I" , but all the following considera- 
tions refer to this picture. 



5. Time Evolution 

To approximate the evolution of the system for a 
short time step, we can truncate the corresponding 
Dyson series 

Mr)) « (l - ^i(r) - ± Mr)) l^(0)> • (17) 

This is a truncation of second order, in which the J/'s 
are the time ordered integrals that occur in the Dyson 
series [8] 

Ui(t)= fdr'VV), (18) 
Jo 

t> 2 (r) = At'V{t) f At"V{t") . 
Jo Jo 

According to the hermiticity of V(t), Ui{t) should be 
hermitian too, which is not the case for £^2(1"). Ui(r) 
has the same off-diagonal form as V(t) whereas U2(t) 
has here a block diagonal form according to the inter- 
action matrix. (To further simplify notation we do not 
write the r dependence of the U's explicitly. Further- 
more, we omit the time dependence of the wave func- 
tion, if it refers to the initial state, i.e., \tp(0)) ■= \ip))- 
As mentioned above we are interested in the time 
evolution of the probability to find the system in its 
excited state W ex (r), or ground state W gT (r), respec- 
tively. Initially we consider W /cx (r). Neglecting all 
terms of higher than second order (products of Ui and 
U2 as well as terms proportional to £/f ) and taking the 
special off-diagonal block form of the interaction into 
account, we get from (17) and (13) 



w c » = + { <<n w x > (19) 
+ ± <v gr it>i 2 iv> sr > - ^ <v> ox i(t>2 + mwn 

W gI (r) = ^ sr |^ sr ) + j t WWA^) (20) 

+ ± <<rwiV x > - ^ {r'm + uiwi 

The strict overall probability conservation requires 
W cx (T) + W gI (r) 

= (V ex (r)|V ex (r)> + (V gr (r)|V 8r (r)) 

= <<nv> ox > + <v> gr i^ gr > = 1 • (21) 

Since the normalization is already fulfilled in the zero 
order, all higher orders must vanish. Obviously the first 
order vanishes automatically. Thus, exploiting (21), for 
the second order of the sum of (19) and (20) we find 

W\(u 2 + ulwi = <vHt>i IVH , (22) 

{^\{u 2 +ut)\^) = {^\m\n ■ 

Inserting this into (19) and (20) yields 
^ cx (r) = (^ x |^ cx ) + ^^ gr |f>i|^ x ) 

-i<VHW sr > 

+ ^W> Sr |^IV> Sr > - ^ cx |t>i 2 |</> OX > (23) 
W^{r) = {^\^)+'^\Ui\^) 
- ^ gr |^ CX > 

+ ^<in#?iin - ^\m\n . (24) 

For an exact evaluation of the right hand side one 
would need to know the |V>° X ), \ip sr ) in detail. But 
rather than doing so, we replace now the actual quan- 
tities by their corresponding Hilbert space averages, 
according to the approximation scheme explained in 
sect. 3. The appropriate set of states over which the 
average has to be taken here is the set of all states 
featuring the same excitation probability (^ ex |V> ox ). 
Some justification for this replacement has been given 
in sect. 3, some comes from the numerical results, sect. 
9, and for full detailed justification, see [9,10]. Since 
not only the justification of the replacement but also 
the actual computation of the respective are beyond 
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the scope of this text we again refer the interested 
reader to [9,10]. Here, we only want to give and discuss 
the results: 



r |^i|V> ex >] = 



1(7x1^)1=0 



(25) 



'■ ST \u 2 i\r r n 



where 



denotes the Hilbert space average and 



Tr ex ( gr ){. . . } the trace over the upper (lower) subspace 
of the operator. 

Plugging those results into (23), we get 



VT x (r) 



(26) 



h 2 Ni 



W gI (r) = 



h 2 N§ 

(0) 



Tr gI {U?} 



v n'Ni 

Now we have to analyze those traces in more detail. 
We will do this explicitly for the upper subspace but by 
simply exchanging the indices, the result will be valid 
for the lower subspace as well 

N c N c 

Tr cx {i>i} = j^(j\m\j) = E|&iJ>f • ( 2? ) 

3 = 1 3=1 

Here j runs over the eigenstates of Ho in the upper 
subspace (note this corresponds to the lower "band" of 
the environment). The object that is summed over here 
is evaluated in the literature in the context of Fermi's 
Golden Rule, 



lrn-\| 2 VM/-it>i ■J 2 4ain2 (l t *.iT) 

\ Ul b)\ = <*l^b> -jr 



l 



(jAE c - iAEi) 



(28) 



(29) 



with 
see Fig. 3(a). 
6. The Linear Regime 



Our arguments, including the conditions we have 
to impose on the model, follow now closely the ones 
brought forth in the context of Fermi's Golden Rule. 




Fig. 2. The function /(w) defined in (31). 

The summation in (28) consists of two different 
terms: the transition elements of the interaction ma- 
trix and a weight /(w). The spacing of different w;j is 
given by 

AEf _ Se 



Aw = uiij — Wj+i : j = 



(30) 



where we have used (10). The function 

/H 



sin 2 (iwr) 



(31) 



is basically a peak at w = 0, with the width Slj — 
4it/t and a height of /(0) = r 2 /4. The area under the 
function /(w) is A = tit/2 (see Fig. 2). This means 
the peak gets higher and narrower as r increases (see 
Fig. 3). 

The height of the peak grows with the square of the 
time r , the area under / only linearly with r. One could 
thus expect two different behaviors: the square- and the 
linear regime. At the very beginning, the peak is very 
broad and therefore much broader than the "band" 
width Se divided by h. In this case we expect that the 
sum grows with the square of r, because all terms are 
near the maximum of the peak (see Fig. 3(b)). We 
choose some n such that the width <5w(tl) of /(w) has 
approximately the same value as the "band" width Se 
divided by h 



<5w(n) 



(32) 



47T Se 4nh 
n n Se 
The terms are distributed over the whole width of the 
peak and we expect that the sum grows proportional 
to the area under the peak, thus linearly in r (see 
Fig. 3(c)) . In this case and if, furthermore, the func- 
tion / does not change much over many summation 
steps Aw, i.e., if 



a Se . . . Se 

Aw = — « M n) = T 



Nt > 1 (33) 
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Fig. 3. Summation of transitions in (28): (a) Matrix el- 
ements to be summed up. (b) r < Ti: almost all terms 
are around the maximum of the peak (square regime) . (c) 
r Ri n : the terms are distributed over the whole peak (lin- 
ear regime), (d) t S! T2: only a few terms arc within the 
peak (break-down of the approximation). 



the summation averages out the different elements of 
the V-matrix in (28). Therefore the sum may be ap- 
proximated by the average of the interaction matrix 
element A§ times the integral over f(u) according to 
cd. The average of the interaction matrix element is 



= 1 2 = 1 



(34) 



For (28) we then get 



Ui\j) 



Auj y ' Auj 



(35) 

where we have used that the area under f(ui) is A = 
7rr/2, as mentioned before. 

The approximation done so far breaks down later 
at some time r 2 , when the peak gets too narrow (see 
Fig. 3(d)), i.e. the width is smaller than the summation 
displacement Alo 

4irHNf 



c / \ 47T Se 
Mr2) = _=Au, = — 



T2 = 



Se 



(36) 

Thus (35) is a valid approximation only for ri < r < 
t 2 , which is the linear regime. 

Hence, plugging (35) into (27) yields 



27rAgfiiV 1 c Af c r 
Se 



(37) 



Since this expression is symmetric under exchange of 
the upper and lower subspace, the corresponding ex- 
pression for the lower subspace reads 



Tr gI {u?}=J^\u 1 \i)\ 



Se 



(38) 



Inserting (37) and (38) into (26) yields 

W cx (r) = VT x (0) + CtN$W si {0) - c*T7vrw ox (o) , 

(39) 

W s \t) = W er (0) + CtN c l W c *(0) - Ct 7V c W/ gr (0) , 
where we have abbreviated 
2ttAo 



Seh 



C . 



(40) 



Equations (39) describe, within the discussed limits, a 
short time step starting from any initial state, not nec- 
essarily an eigenstate of H. Since they directly connect 
the probabilities W c *(0), W SI (0) of the initial state 
with those of the state reached after time r, we can now 
iterate these equations under some specific conditions. 
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7. Conditions on Model Parameters 

Before iterating the above equations (39) , one should 
again check the pre-conditions for the short time step 
equation derived so far. We have only considered terms 
up to second order, and we can only iterate after a time 
step of length n . Thus we have to make sure that the 
considered second order terms are still small compared 
to 1 after n, to justify the dropping of higher order 
terms. Therefore we must check that e.g. 



CtN§\ 



\ 2 



(A£§)2 Ni 



TT7 « 1 



(41) 



where we have used (10). In complete analogy we get 
for the other term of second order 



cvivri 



8^ 



\ 2 



< 1 . 



(42) 



If these two conditions are fulfilled the "linear regime" 
is reached while the truncation to second order is still 
a valid description, and we can iterate (39) after some 
time t > Ti . Obviously the linear regime is reached the 
faster the more levels the environment contains. 

However, if we want to use the above scheme (39) 
we should make sure that we iterate before the linear 
regime is left again, i.e., before T2. Therefore we must 
consider the second order terms at T2 (36) compared 
to one. Note that T2 differs for the two terms of second 
order, in (36) we only argued for one of the two energy 
"bands" in the environment. Thus, the case for which 
iterating (39) is the best description we can possibly 
get is 



W/ cx ((n + l)r)-iy cx (nr) _ 

T 



W gI ((n + l)T)-W el (nr) _ 

T 

CNlW^{nT) - CNZW gI (nr) . 
Or, in the limit of r being extremely small 
dW e 



(45) 



dt 
dW gr 
dt 



= CNoW gI - CNlW™ 
= CN'[W CX -CN^W gI 



(46) 



This evolution equation for the probabilities obviously 
conserves the overall probability. We have obtained a 
rate equation for the probabilities to find the system 
in the upper respectively lower level. 

The solutions of the equations (46) describe simple 
exponential decays, with exactly the same decay rates 
one would have gotten from Fermi's Golden Rule. A 
solution for the considered system being initially en- 
tirely in the exited state reads 



W gI {t) 



Ni + N§ 
Ni 



-C(JVj+JVj)t 



iV c + Ni 



(47) 



iV c + Ni 

The equilibrium values reached after very long times 



W ex (oo) 



N§ 



W g 



Ni 



Ni + iV c v ' iV c + Ni 



(48) 



which are exactly the same as the ones derived in [9] 
for the equilibrium state of a system with an energy 
exchange coupling to a possibly non-Markovian envi- 
ronment. 



CriVo c U T2 = 8^ 2 



CtNI\ 



8tt^ 



(AEiY 



> 1 



> 1 



(43) 
(44) 



8. The Rate Equation 

If the above conditions are fulfilled, iterating (39) 
yields 



9. Numerical Results for the Relaxation Period 

To check the validity of the theory developed in 
the previous Sections, a model of the type depicted 
in Fig. 1, with a Hamiltonian as described in (11) 
has been analyzed numerically by directly solving the 
Schrodinger equation. The interaction matrix V has 
been filled with random Gaussian distributed entries 
such that 

A '' 1, (49) 



(AEi) 
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Fig. 4. Numerical simulation for the relaxation period. The 
predictions from the rate equation get better with increas- 
ing container system, . 



to ensure that (41)-(44) are fulfilled. Different con- 
tainer sizes have been analyzed, corresponding to JVf = 
50, 100, 200, 400, 800 and N§ = §JVf . For all sizes the 
level spacings AEf , AEq have been kept fixed such that 
for increasing container size the band widths increase. 

With those parameter settings the theoretical pre- 
diction for W SI (t) from (47) is the same for all con- 
tainer sizes. The initial state is always chosen to be a 
pure product state {S s (0) = 0), with the gas-system in 
the ground state and the container-system in a random 
superposition of states from the upper band. 

The numerical results are displayed in Fig. 4. The 
solid line is the prediction from theory. Obviously the 
theoretical predictions are not accurate for "few-level" 
container environments. This is due to the fact that 
the replacement of actual quantities by their Hilbert 
space average is a valid approximation for high dimen- 
sional Hilbert spaces only. Furthermore, for the few- 
level cases the iteration step times that have to be 
longer than n are rather long, because already n is 
long. This means that the recursion cannot really be 
replaced by the differential equation in Sect. 7. This es- 
sentially shifts the corresponding curves to later times, 
compared to the theoretical prediction. All those ef- 
fects vanish if the container system becomes sufficiently 
big. The simulation for JVf = 800 is obviously in good 
agreement with the theoretical prediction. 



10. Summary and Conclusion 

We considered a two-level system coupled to an envi- 
ronment specified only by two resonant energy bands. 
We showed that for systems of this type (or rather any 
bi-partite quantum system) the entropy of the consid- 
ered system cannot increase without increasing system- 
environment correlations, if the full system evolution 
is unitary. 

Further more we solved the Schrodinger equation 
for the full system starting with an initial pure prod- 
uct state. Under some conditions concerning interac- 
tion strength, band width and state density of the 
environment-system we find a merely statistical energy 
transfer or relaxation process, that may simply be de- 
scribed by some transition rates. This result has been 
derived on the basis of a theory involving Hilbert space 
averages and confirmed numerically. 

We thus conclude that for statistical decay behav- 
ior of exited states in quantum systems neither ther- 
mal nor infinite baths are necessary, just as well as the 
factorizing condition seems neither tenable nor indis- 
pensable. 
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